Mon. Not. R. Astron. Soc. 000, [l]{7] (2008) Printed 13 September 2009 (MN IATgX style file v2.2) 



Disks in the sky: A reassessment of the WMAP "cold spot 

Ray Zhang 1 * and Dragan Huterer 1 ! 

1 Department of Physics, University of Michigan, 450 Church St, Ann Arbor, MI 48109 



13 September 2009 



ON 
O 
O 

(N 
Oh 



O 

U 

6 

> 

oo 

00 
ON 

cn 

od 
o 

ON 

o 

> 

X 
C3 



ABSTRACT 

We reassess the evidence that WMAP temperature maps contain a statistically sig- 
nificant "cold spot" by repeating the analysis using simple circular top-hat (disk) 
weights, as well as Gaussian weights of varying width. Contrary to previous results 
that used Spherical Mexican Hat Wavelets, we find no significant signal at any scale 
when we compare the coldest spot from our sky to ones from simulated Gaussian 
random, isotropic maps. We trace this apparent discrepancy to the fact that WMAP 
cold spot's temperature profile just happens to favor the particular profile given by 
the wavelet. Since randomly generated maps typically do not exhibit this coincidence, 
we conclude that the original cold spot significance originated at least partly due to 
a fortuitous choice of using a particular basis of weight functions. We also examine 
significance of a more general measure that returns the most significant result among 
several choices of the weighting function, angular scale of the spot, and the statistics 
applied, and again find a null result. 
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INTRODUCTION 

Cosmic microwave background (CMB) maps have been 
studied in detail during the last few years. These stud- 
ies have been motivated by the remarkable full-sky high- 



resolution maps obtained by WMAP (Bennett et ah 2003 
Spergel et al.|[2007 l, and led to a variety of interesting and 



unexpected findings. Notably, various anomalies have been 
claimed pertaining to the alignment of largest modes in the 



CMB (Tegmark et al.| 


2003; 


de Oliveira-Costa et al. 


2004 


Copi et al. 2004, Schwarz et a 


. 2004; 


Land & Magucijo 


2005 


Copi et al. 2006 


Abramo et al.|2006 


1, the missing power on 


large angular scales (Spergel et al. 2003 Copi et al. 2007 



Copi et al.||2008p , and the asymmetries in the distribution 



of power (|Eriksen et al. 2004 Hansen et al. 2004 Bernui 



et al.|[2006 Hajian 2007). In the future, temperature maps 



obtained by the Planck experiment, and large-scale polariza- 
tion information ( Dvorkin et al.|2008[ ) may be key to deter- 
mining the nature of the large-scale anomalies. For a review 
of the anomalies and attempts to explain them, see |Huterer| 
( 2006] ). 

Several years ago, Vielva et al. (20041 reported an 



anomalously cold spot in the WMAP microwave signal: kur- 
tosis of the distribution of spots (defined using Spherical 
Mexican Hat Wavelet weight functions) is unusually large 
on scales of about 5°, at < 0.5% significance. The authors 
also noted that the result is driven by a a cold spot in the 
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southern hemisphere, at (£,b) = (-57°, 209°). The finding 
has been confirmed and further investigated in |Cruz et al.| 
(2006), who found that an equally cold or colder spot of 
this size is expected in < 1% of Gaussian random, isotropic 
skies, as well as Mukherjee & Wang (20041; Cruz et al 



(2005) 
p006 > 



McEwen et al.|(|2005[);|Cay6n et ai. | J2"005F ; 
Rath et al.|(|2007|;|Naselsky et al.|p007) 



et al. (20081; Rossmanith et al. (2009), some of whom also 



Cruz et al. 



Pietrobon 



studied the spot's morphology. The plot further thickened 
when Rudnick et al. (20071 claimed that there is a corre- 



sponding cold spot (underdensity in galaxy counts) — in 
the NVSS radio survey, and at roughly the same location 
as the CMB cold spot; however, this particular claim was 



shown by Smith & Huterer (20081 to be an artifact of the 



a posteriori statistics and the particular way NVSS data 
had been analyzed. Nevertheless, the CMB cold spot re- 
mains a much-studied topic and the source of investigations 
of whether exotic physics could be the cause. 

Perhaps surprisingly, nearly all of the works so far con- 
sidered searches for the cold spot using the same basis func- 
tions — Spherical Mexican Hat Wavelets (though with a few 
exceptions — Pietrobon et al. ( 2008 1 used needlets, while 
Rath et all 12007b and IRossmanith et all (I2009I) used the 



scaling indices). The only variation in the different analyses 
was in the choice of the statistics that was applied to the 
wavelet-based weights. 

Here we set out to check the evidence for the cold spot 
using different, and arguably simpler, set of weight func- 
tions. We reassess evidence for the "cold spot" using circu- 
lar top-hat weights (i.e. disks) of arbitrary radius R. We do 
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so in order to verify findings that relied on wavelets, and 
more generally to investigate the robustness of the signal. 
We also check results using simple Gaussian weights, find- 
ing results consistent with those with the disks. We then 
investigate the source of this apparent discrepancy with all 
of the previous work that used wavelets, and find that the 
cause of the discrepancy is the specific temperature profile 
of the cold spot which just happens to favor the profile of the 
Spherical Mexican Hat Wavelet. In addition to the choice of 
the spots' weight function, the original claim refers to the 
angular size of the spot of ~ 5° that is also chosen a poste- 
riori. We investigate the effect of these choices by defining 
a "superstatistic" measure that combines several previously 
considered statistical measures of coldness and the associ- 
ated choices of the spot size and weight functions, and find 
that the claimed spot (or any other spot in our sky) is not 
unusually significant using this new measure. 



STATISTICS AND MAPS 
Weight statistics 

The top-hat weights are familiar from structure formation 
(where they are used in the definition of the amplitude of 
mass fluctuations over some scale R, for example) and effec- 
tively represent another statistic to study the cold spot. We 
define the disk top-hat weight of radius R as 

D(r)=A dis w(R)[G(r)-G{r-R)], (1) 

where Q(x) is a Heaviside step function and Aiisk(-R) = 
(2tt(1 — cos(_R))) -1 / 2 is defined so that 



D(rfdQ. = 1. (2) 



Note however that the normalization yldisk(-R) is unimpor- 
tant for finding the coldest spot since we only do relative 
comparisons of temperatures in disks on the sky. The top 
hat-weighted temperature coefficients are given by 



T disk (r;R) = I dQ'T(f')D(a;R), 



(3) 



where f = (8, (j>) is the location of a given spot, f ' = (6' , <f>') 
is the dummy location on the sky whose temperature we 
integrate over, and a — arccos(f • f') is the angle between 
the two directions. 

The Gaussian weights that we use are defined equiva- 
lently. The weight functions are 

„2 



G(r) = .4 GausB (i?) exp -4 In 2 



R< 



(4) 



so that the full width at half maximum of the distribution 
is equal to R. The weighted temperatures are given by 



T GauBS (f;R) = J dQ'T(r')G(a;R). 



(5) 



Finally, the corresponding procedure applied to the 
wavelets is as follows ( |Cayon et al.| 2001 Ma rtinez-Gonzaiez| 
et al.|2 002 ) . The Spherical Mexican Hat wavelets are defined 
as 



<f{9;R) = A mv (R) 1 + 



2 - 



R 



exp 



\2R 2 J' 
(6) 



where y = 2tan(6 l /2) and 



2nR A 



R 2 R 4 



-1/2 



(7) 



so that J dQ^ 2 (6; R) — 1 over the whole sky. We can now 
define the continuous wavelet transform stereographically 
projected over the sphere with respect to ty(9;R), with T 
being the CMB temperature: 



Twav(z;i?) = j dQ'T{x + p')^(e'-R), 



(8) 



where x — » (8,<f)) and pJ — > (9',<j)') are the stereographic 
projections to the sphere of the center of the spot and the 
dummy location, respectively, and are given by 



2 tan - (cos <f>, sin < 



2 tan — (cos cf> ', sin </>') 



(9) 



(10) 



Martinez-Gonzalez et al. (20021 for details. To work in 



terms of purely spherical coordinates, we center the spot 
location to the north pole of the sphere, and rewrite the 
above as 



(f;R) = / dn'T(f')^(a;R), 



(11) 



where M(f') is the mask, defined to be 1 for pixels within 
the mask and and for those outside of it. As the wavelet 
is effectively zero for a values greater than ~ 4 times the 
radius, we can carry out the integral by using the Healpix 
command query _disc to find all pixels within a circle of that 
radius from the wavelet center. 

To account for the masked parts of the sky, at each spot 
location f we first calculate the "occupancy fraction" 



N{r;R) = J dQ' 'M(r')* 2 \0; R). 



(12) 



We only include results for spot locations f for which 
N(f;R) > 0.95. Additionally, we do not include individ- 
ual pixels that have M(f') < 0.9 in order to limit biases 
due to masking (partially masked pixels come about after 
degrading maps to a lower resolution). As discussed further 
below, we tested our procedures by using a higher occupancy 
fraction and found consistent results. 



Maps 



We use WMAP's five year maps in our analysis (Hinshaw 



|et al.|2 009). Following |Vielva et al] ( |2004[ ), the fiducial map 
we use is the coadded foreground-cleaned map 



T : 



J2ll 3 T r (l)W r (i) 



(13) 



where T is the coadded temperature, determined from the 
weighted sum of temperatures T r of each individual radiome- 
ter r e {Ql, Q2, VI, V2, Wl, W2, W3, W4}, divided by the 
total weight. The weights at each pixel for each radiometer 
are w r (i) = N r (i)/a 2 , where N r (i) are the number of effec- 
tive observations at the pixel, and oy is the noise dispersion 
for the given receiver. 

This coadding was performed on maps at resolution of 
A s idc = 512 (~ 8'), then the KQ75 mask was applied. As 
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Figure 1. Significance in "sigmas" of the cold spot temperature S coldest (.R) (left panel) and kurtosis of the distribution of S(R) (right 
panel) for the three choices of weights that we have examined: disks (black), Gaussians (blue) and wavelets (yellow). Here "lcr" means 
< 32% likely, "2cr" means < 5% likely, and "3<x" corresponds to < 0.3% . Note that the scale R, shown on the x-axis, is defined separately 
for each choice of the weighting function so that the corresponding fair comparison can be made; see text for details. The results show 
that only the wavelet-based temperature cold spot and kurtosis deviate significantly from the Gaussian random expectation, and that 
the disk or Gaussian ones do not. 



mentioned earlier, spots with more than 5% of the weighted 
area masked (N(r;R) > 0.95) were not used. 

The locations of centroids of spots are chosen to be cen- 
ters of pixels in N B ide = 32 resolution; therefore, we examine 
iVpix = 12 N? dc ~ 12, 000 spots on the sky. In order to calcu- 
late the spots' weighted temperatures, however, we analyze 
the coadded map at the N s i dc = 128 (~ 0.5°) resolution, 
which is sufficiently high to lead to converged results for 
R ;> 2° spots, yet sufficiently low to be numerically feasible. 

The results of our analysis were then compared to 
10,000 randomly generated Gaussian full sky maps, with 
the same methodology applied. The skies have been gener- 
ated using the Healpix facility synf ast, and used as input 
the power spectrum determined in the WMAP 5-year anal- 
ysis ( Nolta et al.||2009 |. The maps were then smoothed by 
a Gaussian with FWHM = 1° to match the WMAP proce- 
dure. 



Significance statistics 

The principal statistic that we use is the temperature of 
the coldest spot divided by the standard deviation of the 
distribution of all spots 



S diak (f;R) 



TIC 



(r,R) 



0"disk(-R) 



(14) 



and equivalently for the Gaussian weights and the wavelets. 
Here <Jdisk(R) is the standard deviation of the distribution 
of all spots in a given map, while T d °^ eBt (R) is the coldest 
spot in the distribution. Note that the distribution of spot 
temperatures is not Gaussian as we noted earlier, but this 
is irrelevant for us; we scale T by a in Eq. ( |14[ ) in order to 
account for small (~ 10%) differences in the overall level of 
power in spots of characteristic size R in the different maps 
— in effect, o"disk(-R) provides units in which to best report 
the coldest temperature. 

Computing the significance of our statistic Sdisk(r;R) 



is then in principle straightforward: we compare it to values 
obtained from simulated Gaussian random maps and rank- 
order it; the rank gives the probability. 

In addition to the cold spot significance, we follow 



Vielva et al. (2004) and Cruz et al. (2005) and consider 



the kurtosis of spots in a given map. The kurtosis is simply 
related to the fourth moment of the distribution of the spots 



K diBk (R) 



1 T£r m T iUk (f i ,R) 4 



pots 



CTdisk(ii) 4 



(15) 



and equivalently for the Gaussian weights and the wavelets. 



RESULTS 

Wavelet weighted spot 

We first make sure that we reproduce the cold spot results 



of Vielva et al. (20041 and Cruz et al. (20051. For Spheri- 



cal Mexican Hat Wavelets with R = 5°, we find the center 
of the coldest spot in the five-year combined cleaned map, 
is at coordinates (t, 6) = (—57.7°, 209.3°) (corresponding 
to spherical coordinates (9,0) = (147.7°, 209.3°)). In gen- 



eral agreement with the Cruz et al. (20071 results, we find 



that only (0.99 ±0.10)% of simulated statistically isotropic, 
Gaussian random maps exhibit a more significant cold spot 
(i.e. a more negative value of S^$ eBt (R)) for this value of 
R. Here and throughout, the error bars account for the fi- 
nite number of (N — 10, 000) simulated maps; we quote the 
standard margin of error which, for a fraction p of a total 
of N events, is given as a(p) = \/p(l — p)/N. Moreover, we 
confirm that while the variance and skewness of the distri- 
bution of S^ CBt (R) from synthetic maps, the kurtosis at 
R = 5° is high at the (1.44 ± 0.12)% confidence. 

Reporting the significance only for the R = 5° may be 
unfair, however. To address this, we consider the range 2° sC 
R 5j 8° in steps of 0.5° (the lower bound is set to correspond 
to spots significantly larger than smoothing of the maps of 
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Figure 2. Left panel: disk and wavelet weights, together with the azimuthally-averaged CMB temperature profile as a function of radial 
distance from the (wavelet-found) center of the cold spot. Right panel: contribution to the weighted temperature using disks and wavelets; 
the blue (red) shaded areas show the difference by which wavelets (disks) dominate in the given ranges of distance, while the dashed line 
shows the cumulative difference between the two. Note that all lines in both panels have arbitrary normalization, although the disk and 
wavelet lines are consistently compared using their fiducial normalizations from the text. 



roughly 1°). The significances of these results are given as 
black lines in the two panels of Fig. [J] The significance of 
the wavelet-determined cold spot peaks around R = 5-6°, 
while the kurtosis is significant in the range 3° ^ R ^ 5°. 



Disk and Gaussian weighted spot 

We now repeat the same analysis with the circular weights. 
First we confirm that the coldest disk-weighted spot in the 
co-added smoothed Q-V-W map is at nearly the same loca- 
tion as the wavelet-weighted spot, at (I, b) = (—57.4, 208.0). 
The result is similar in the foreground cleaned Q-V-W map. 
However, the disk-weighted spot is not unusually cold: at all 
scales R between between 2° and 8°, the statistic S^°l^ eBt (R) 
is not unusually low relative to expectations from Gaussian 
random maps; see the left panel of Fig. [l] We find the same 
results for the kurtosis of the distribution of Sdisk(r;R) — 
as shown in the right panel of Fig. [I] these distributions fall 
well within the expectation on all scales we examined. 

Surprised by these results, we have repeated the same 
tests with Gaussian weighting, where the full-width half- 
maximum (FWHM) of the Gaussian weight has been set 
equal to the scale R; this way we ensure that a large fraction 
(about 94%) of the weight is applied within the radius R. 
The results are similar as for the disks in that they are not 
significant; see again Fig. [I] In particular, the coldest spot 
is most significant at R = 4-5°, but even there only at the 
1-cr (25%) level, while the kurtosis is not significantly large 
or small at any scale. 



1 These tests are computing-intensive, and we were forced to 
compare to simulated maps at A^j^e = 32 (rather than N 3 ^ e = 
128) for wavelets for 6.5° SC R ^ 8°; we have checked at lower 
values of R that the results at the two resolutions are in good 
agreement. 



The density profile of the cold spot 

The question is obvious: Why was the cold spot so signifi- 
cant for wavelets, but not so much for disks? To address this, 
we show the disk and wavelet weights, together with the 
azimuthally-averaged CMB temperature profile, as a func- 
tion of radial distance from the (wavelet-found) center of 
the cold spot in the left panel of Fig. [2] In the right panel, 
we show contribution to the weighted temperature using 
disks and wavelets, as well as the cumulative difference be- 
tween the two. This figure shows the case of R = 5° which 
approximately maximizes significance of the wavelet-based 
cold spot. 

The azimuthally averaged density profile of the temper- 
ature is about equally distributed between zero and 5° (that 
is, the blue area in the right panel is about the same as the 
red one up to 5°). However, at distance beyond the edge 
of the disk of 5°, the wavelet accumulates more weight as 
seen in the right panel. The reason is shown with the curve 
labeled "averaged CMB temperature" in the left panel: the 
CMB profile goes from negative to positive with increas- 
ing radius from the center of the (wavelet-based) cold spot, 
precisely favoring the wavelet profile that has roughly the 
opposite behavior. 

What is the likelihood of this conspiracy that the tem- 
perature profile of the coldest spot mimics the shape of 
the wavelet? Using Gaussian random maps, we estimate 
the likelihood that a given map has the wavelet-determined 
cold spot is more significant than the disk-determined cold 
spot by at least as much as in WMAP where \S^ est (R) - 
5'J i ° 1 c dcst | = 1.5. We find that the wavelets are more significant 
that the disks by at least this margin in only (1.89±0.13)% of 
the random maps (while the disks are as or more significant 
in only (1.96 ± 0.14)% cases). From this, we conclude that 
typical Gaussian random CMB maps do not show increased 
significance of the wavelet determined cold spot, relative to 
the disk-determined one, to the same extent as our sky does. 
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Cold spot statistic 


(R = 5°) 


Kurtosis statistic (R 


= 5°) 




Wavelet 


Disk 


GetUSSicLH 


Wavelet Disk 


GcLiissicin 


Value of statistic 


-3.21 


-4.54 


-3.52 


0.58 -0.23 


-0.06 


Significance in % 


0.99 ±0.10 


46.3 ±0.5 


25.1 ± 0.4 


1.44 ± 0.12 69.3 ±0.5 


46.5 ± 0.5 



Table 1. Statistics and its significance for coldest spot and kurtosis of spots, evaluated on scale R = 5°, for the three weights we 
considered. For the coldest spot, the statistics is S, defined in Eq. l |14| l, while for kurtosis, the statistic is just its value, defined in 
Eq. jl5}. The error bars are Poisson and reflect the finite number of the simulated Gaussian random maps to which we compared 
WMAP. The results are robust to variation in the choice of the WMAP map, or radius R, as discussed in the text. 



Robustness 

We have tried varying a number of details, with the following 
results: 

• In addition to the coadded foreground cleaned Q±V±W 
map, we also used the coadded Q±V+W map, the coadded 
V±W map, and the coadded foreground cleaned V±W maps 
available from WMAP. All of the significances are very close 
to the foreground QVW map value; we have checked that 
smaller scales (1-1.5°), which were not in the range we pre- 
sented in the final analysis, would be somewhat discrepant; 
this is not surprising given that the maps are smoothed to 
1°. 

• To test the effects of the resolution of the map, we vary 
the resolution from N B id e ~ 128 (~ 0.5° pixels) to N B id c = 32 
(~ 2° pixels), with 16 times fewer pixels. We again find con- 
sistent results except at small scales, R < 2°, which makes 
sense since pixelization is expected to play a role only when 
pixel size becomes comparable to the spot scale R. 

• To ensure that our finite step size has not accidentally 
overlooked a cold spot, we refine the resolution in our search 
for the coldest spot in WMAP by querying at every pixel in 
an A'side = 512 map within 3° of the center of the reported 
cold spot; this stepping size is effectively ~ 256 times higher 
than before. While an increase in the temperature of the cold 
spot is entirely expected, we find that this increase is small 
enough not to appreciably change the significance results for 
all three choices of the weight function. 

• To test our prescription for dealing with partially 
masked spots, where we only analyze spots that have the 
"occupancy fraction" N(f;R) > 0.95 (see Eq. ( |12[ )), we re- 
peat the analysis with the minimum occupancy number of 
0.98. While the resulting number of spots retained in the 
analysis is now much smaller, decreasing by between tens of 
percent (for spots at R — 2°) to about a factor of 10 (for 
R = 8° disks), we find results generally consistent with our 
fiducial case: the statistic S and kurtosis calculated using 
the wavelet weights are significant, while the same statistics 
calculated using the disk and Gaussian weights are not. 



More general tests 

It is clear that a posteriori choices were made in the original 
claims for the existence of the cold spot — in addition to 
the choice of the weighting function (which is the principal 
subject of this paper), the moment of the distribution of 



spots (kurtosis) and spot scale (5°) have been called out 
after noticing that they are unusuarl In contrast, variance 
and skewness of the spot distribution, or kurtosis and scales 
larger or smaller than ~ 5°, do not show departures from 
expectations based on Gaussian random isotropic maps, as 
we have checked as well. 

We have investigated how results change with more gen- 
eral tests as follows. We have formed a "superstatistic" de- 
fined as maximum significance of either variance, skewness, 
kurtosis, or coldness (the last two being defined earlier in 
this section) over any scale R or weight function set W 

(16) 



max {P(Stativ(7?))} 

R, W, Stat 



wher^] 

R e 
W G 
Stat G 



{2°,2.5°,...,8°} 
{wavelet, disk, Gaussian} 
{Variance, Skewness, Kurtosis, 5*} 



(17) 



and where each probability P was individually calculated 
relative to Gaussian random isotropic maps as described 
earlier. [The S and kurtosis statistics have been defined in 



Eqs. ( 14 1 and ( 15 I, while the variance and skewness are de- 
fined analogously to kurtosis.] Note that we define P to cap- 
ture the possibility that the statistics in question is either 
small or large relative to expectation; in other words, we 
adopt the minimum of r and (100% — r), where r is the rank 
of the statistics relative to simulated maps. 

We find that the value of the iS sup cr statistic for the 
WMAP cleaned QVW map is 0.54%, and this value is at- 
tained by the kurtosis statistic Swavcict at scale R = 3.5°. 

2 While observation of some of the other anomalies mentioned in 
the introduction was technically also a posteriori, those anomalies 
had to do with special scales (e.g. largest observable scales on the 
sky) or directions (e.g. the ecliptic, which the telescope pointings 
preferentially avoid). In contrast, there appears to be nothing 
special about kurtosis of the spot distribution, or scales of ~ 5°. 

3 For wavelets with 6.5° R $C 8° the calculations at N s i^ c = 
128 were unfeasible due to the large number of pixels to keep 
track of. While we have checked that N a i^ e = 32 wavelet results 
are similar to N s i^ B = 128 at R < 6.5°, and show no significant 
results at larger scales, for consistency we decided to quote the 
superstatistic results at N s ^ c = 128 and consider the 2° Sj R ^ 
6.5° scales for the wavelets, and 2° $J R $J 8° for the disks and 
Gaussians. 
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Superstatistic value (%) 

Figure 3. Significance of the superstatistic iS sup or, denned in 
Eq. |16) , as a function of its value. The line shows results cal- 
ibrated from our simulated maps, while the point denotes the 
value (and the unremarkable significance) from WMAP. 



is a chance event, since only 5% of the Gaussian random, 
isotropic skies exhibit equal or more significant discrepancy 
in favor of the wavelets. Moreover, we found that the result 
is insensitive to the choice of the map or the statistic used 
for the cold spot. 

Motivated by these findings, we also examined signif- 
icance of a more general measure - which we called the 
"superstatistic" - that combines the various choices of the 
weighting function, spot size, and statistics, and returns the 
most significant choice consistently for each map. We again 
find a null result; the WMAP superstatistic is low only at 
~ 20% level relative to Gaussian random and isotropic sim- 
ulated maps. 

Therefore, we find no compelling evidence for the 
anomalously cold spot in WMAP at scales between 2 and 
8 degrees. The existing evidence apparently hinges on the 
particular choice of the weight functions to define the spot 
(Spherical Mexican Hat Wavelets) and their scale [R ~ 5°). 
While our conclusion may sound like a depressing null re- 
sult, we are upbeat about future tests with WMAP (and 
soon, Planck) to uncover and test unexpected features and 
anomalies. 



However, this value is not too unusual: we find that 23% of 
the Gaussian random skies have a smaller value of 5 sup er; 
see Fig. [3] 

To test the robustness of this result, we consider an 
alternative, more restricted, definition of the superstatistic 
where only the wavelet weights are considered, but where 
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we still varying scale R and statistic Stat; see Eq. ( 18 1. Here 



we effectively assume that, for whatever reason, wavelets are 
the preferred weight functions to be used, but we still seek 
to avoid the a posteriori choices of the scale and statistics. 
The new superstatistic is again not statistically significant; 
we find that 15% of Gaussian random skies show greater 
significance. 

Thus, superstatistic results confirm our earlier conclu- 
sion that less a posteriori tests do not indicate a statistically 
significant cold spot in the WMAP data. 



CONCLUSIONS 

The "cold spot", together with low power at large angles, 
multipole alignments, north-south power asymmetry, has 
been one of the most studied anomalies in WMAP CMB 
temperature maps. So far there have been no compelling 
proposals, cosmological or systematic, that would explain 
existence of the claimed cold spot, which is perhaps not sur- 
prising given that neither its radius (~ 5°) nor its direction 
in the sky are particularly special. 

In this paper, we have investigated evidence for the cold 
spot. While we confirmed its high statistical significance us- 
ing the wavelet basis of weight functions, we did not confirm 
the existence using the disk top-hat, or Gaussian, weights. 
The cold spot is indeed at the same location in WMAP maps 
with the latter two bases, but it is not significant when com- 
pared to expectation based on Gaussian random, isotropic 
skies. 

We traced the apparent inconsistency to the fact that 
the radial temperature profile around the cold spot center 
is such that it favors the wavelet profile; see Fig. [2] This 
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